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A new method to calculate the spectrum using cascaded 
open systems and master equations is presented. The method 
uses two state analyzer atoms which are coupled to the sys- 
tem of interest, whose spectrum of radiation is read from the 
excitation of these analyzer atoms. The ordinary definitions 
of a spectrum uses two-time averages and Fourier-transforms. 
The present method uses only one-time averages. The method 
can be used to calculate time dependent as well as stationary 
spectra. 

PACS numbers: 42.50.Kb, 32. 90. -ha, 32.70.Jz 



I. INTRODUCTION 

Much of the information obtained in laser experiments 
comes in the form of spectral data. In steady state, these 
are supposed to give the Fourier transform of the energy 
level structure of the system under investigation. In this 
manner we have been able to learn about the quantum 
configurations of atoms, molecules and solid state sys- 
tems. 

Today, however, much work is done with pulsed laser 
sources, where the pulse duration samples the evolu- 
tion of the system over times ranging from nanoseconds 
to femtoseconds. Spectral data are still recorded, but 
their significance and physical information is no longer 
straightforward. Thus one needs to reconsider the def- 
inition and computation with due consideration of the 
physical conditions under which the data are obtained. 
Several attempts have been made to modify the Fourier- 
transform definition, the Wiener-Khintchine spectrum 
[Q-H, but here we choose to consider the physical spec- 
trum introduced by Eberly and Wodkiewicz [Q, which 
attempts to emulate the spectral measurements using e.g 
a Fabry-Perot filter in front of the detector. They have 
applied it to the well known case of the fluorescence Mol- 
low spectrum A very similar approach has allowed 
Kowalczyk et al. to consider the fluorescent spectra 
deriving from a laser-excited molecular wave packet. A 
model calculation of such a process has been presented 
by Vinogradov and Janszky |^ . Another approach to the 
molecular situation is presented by Lee et al. 

In this paper we consider the problem of obtaining 
spectral data in an evolving system. Thus only the ac- 
cumulated information is available; the future is still un- 
known. This rules out the use of a full Fourier-transform, 
and only physically manipulated collected data can be 
used. The spectral measurements require a filtering. 



which smeares the signal in time. The frequency-time 
resolution has to obey an uncertainty relation. In ad- 
dition, the exact nature of the transfer of spectral in- 
formation from the system investigated to the detector 
imposes its own limitations. There can be no unique 
spectral definition for time-dependent systems, but we 
can require that all definitions agree when infinite mea- 
surement times are available. 

In a laboratory experiment, the radiation emitted from 
a driven system reaches the detector through a technical 
setup which, in addition to its function as a filter, will im- 
pose its own noise limits. The filtering action provides a 
noisy channel. This is, however, just the situation which 
is described by the term "Open system" We want 

to consider an open systems approach to time-dependent 
spectra. 

A natural frequency-selective detector is a two-level 
atom. This will respond to radiation only within a band- 
width given by its natural linewidth. If this is small 
enough, the spectrum if incoming radiation is resolved 
with this accuracy; this mirrors the quantum fluctua- 
tions of the detector-atom decay. We also let the radia- 
tion reach the detector atoms through a noisy channel, a 
reservoir which can be eliminated. Thus we present a new 
approach to time dependent spectra, which does not seem 
to be directly related to the ones used earlier. The aim of 
the present paper is to introduce this model and compare 
it with the physical spectrum of Eberly and Wodkiewicz 
and the Wiener-Khintchine spectrum for steady state. 

In Sec. II. A we introduce the system we want to inves- 
tigate the spectrum of. We choose a simple one consist- 
ing of three levels only. In fact, we would like to consider 
some real situation like an excited molecule but the com- 
putational burden imposed by the theoretical methods 
is so large, that we have been forced to carry out our 
comparisons on this simple model. 

Sec. II.B.l explains how the stationary Wiener- Khint- 
chine spectrum can be calculated in the Schrodinger pic- 
ture using the quantum regression theorem, this is used 
as a reference in the following. In Sec.II.C.l it is shown 
how the quantum regression theorem can be used to ob- 
tain the physical spectrum of Eberly and Wodkiewicz. 
We will use these results to judge the reliability of the 
method presented in this paper. Sec. II. C. 2 presents the 
theory of our own approach, which relies heavily on the 
work by Carmichael . The derivation leads to a mas- 
ter equation, which is used in a quantum simulation to 
obtain spectra as explained above. The results are pre- 
sented in Sec. Ill and compared with the physical spectra 
in Sec. III.C. Finally Sec. IV contains the conclusions. 
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II. MEASURING THE TIME-DEPENDENT 
SPECTRUM 

A. The system under investigation 

The model system we choose for our investigation is the 
3- level atom shown in Fig. 1. The three levels are taken 
to have the (dimensionless) energies: toi = 0, uj2 = 4.0 
and L03 = 8.0. The levels |2) and |3) decay to level |1) 
by the rates Fi = r2 = 0.1. These parameters are kept 
constant all through the present paper. The level pairs 
1 2) < - > |1) and |3) < - > |2) are coupled through 
lasers with the Rabi frequencies Qi and Q2 respectevely. 
These coupling strengths are varied in the calculations. 
The fact, that the dipole approximation docs not allow all 
transitions considered, lacks significance for the features 
we are investigating in this paper. In a molecular system 
parity may not be a good quantum number, in an atom 
or quantum dot system, the decay F2 may derive from 
higher multipole transitions. 



iJ = ^a;,|z)(zH-^(|l)(2|e--* + h.c) + 

i=l 

^(|2)(3|e--*+h.c). (3) 



B. The measurement scheme and the stationary 
spectrum 

The usual setup for spectrum measurements is pre- 
sented in Fig. Ij. Fluorescence radiation from the three 
level atom is guided to the filter which allows radiation 
with only a certain frequency ujd to go through it. Be- 
hind the filter, a photodetector detects the intensity of 
the radiation coming through. The spectrum can be 
measured by scanning ujjj over an appropriate frequency 
range and recording the relative intensities. The same 
measurement configuration can be used to measure both 
stationary and time-dependent spectra. 



CO, 




FIG. 1. The three level atom used in the calculations. The 
atom has three energy levels oji = 0.0, uj2 = 4.0 and uj^ = 8.0, 
two internal decays Fi = r2 = 0.1 and two lasers Qi and $12- 

Mathematically the dynamics of this kind of system 
can be studied using master equations. We have included 
two different spontaneous decays to the reservoirs so the 
master equation has two decay terms of the Lindblad 
form. The reservoir is taken to be at zero temperature. 
Assuming the Born-Markov and the Rotating Wave Ap- 
proximations (RWA) we find the master equation 

d' ^ II 

J = ~i[H,g] + J2^.,{UQLl - -LlUg--gLlL,) (I) 

1=1 

where Li = |I)(2| and L2 = |I)(3| are Lindblad opera- 
tors. Everywhere in this paper the units have been chosen 
in such a way that h — 1. Lasers are taken into account 
by adding resonant driving terms to the Hamiltonian; the 
laser frequencies are chosen such that 

OJI2 = - ^2- (2) 

In the rotating-wave approximation the Hamiltonian 
becomes 



Fluorescence radiation 




FIG. 2. The measurement scheme. The system is driven 
by a laser. Part of the fluorescence radiation passes the filter 
to the detector. 



1. Stationary spectrum and the Quantum Regression 
Theorem 

Mathematically the spectrum for a stationary system, 
the Wiener-Khintchine spectrum, is defined as 

00 

Swk{oj)^ j dT{a\t)d{t + r))e'^\ (4) 

—00 

where a and are field operators. The definition con- 
tains two-time averages of the field modes. These expec- 
tation values are averages over the degrees of freedom of 
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the system, so in practise the spectrum is calculated from 
two-time system expectation values. 

In the calculation of the spectrum, two-time averages 
of system operators in the Heisenberg picture are needed. 
Simulations are, however, done in the Schrodinger pic- 
ture. The method which allows us to calculate multitime 
averages in the Heisenberg picture using a master equa- 
tion is called the Quantum Regression Theorem (QRT). 

The calculation of the expectation value {02{t + 
T)Oi{t)), T > can be carried out with the QRT us- 
ing the following steps: 

1. Evolve the density matrix g using the master equa- 
tion from the initial time to to the time t. 

2. Form the matrix Cij{t, r = 0) 



Cy(i,T -0)-Tr[e(i) 10 (j|Oi 



(5) 



where Oi is the operator in the Schrodinger picture. 
3. Evolve Cij (t, r) using equation (|^) with the initial 

" ■ 0) 



matrix Cij (t, r 



dC\,it,T) 
dr 



GijkiCki{t, r); 



(6) 



here repeated indices are summed over. Equation (|6|) is 
the master equation (^) in component form. The coeffi- 
cients Gijki are thus determined by the master equation. 

4. Form {02{t + T)Oi{t)) using the matrix elements of 
O2 in the Schrodinger picture ie. 



{d2{t + T)di{t)) ^^2,^JaJ{t,T). 



(7) 



This procedure gives the expectation values for such time 
pairs (t, t + t) where the operator O2 has a larger time 
value than the operator Oi, ie. r > 0. In order to get 
time values in the reverse order, we have to take a com- 
plex conjugate. It is also possible to calculate expectation 
values of the type (di(i)02(i + T)03(t)). The difference 
is that at step 2 we have to calculate 

C,,it,T = 0)=Tr[\t){j\digit)d3] (8) 

and then continue as in the previous case. 



2. The Stationary spectrum for the three level system 

Without performing any calculations, we expect the 
following kind of structure of the stationary spectrum of 
spontaneous emission for the three level system. When 
the laser amplitudes are small, the spectrum should have 
two peaks whose relative intensities depend on the pa- 
rameters Fi and F2. If Fi is large, then only level two 
should decay significantly and the peak from transition 
from level three to one should be small. When the laser 
intensities are increased, we should see some kind on Rabi 
splitting in the spectrum, ie. the two peaks are expected 
to display substructure. 



The detection operator used in the calculations is 
= Fi|l)(2|+F2|l)(3|. 



(9) 



The two-time average {O'' {t)0{t+T)) has been calculated 
and the spectrum is obtained by integrating this expec- 
tation values with respect to r. In Figs. || and ^ the 
stationary spectrum of our three level system is shown 
for two different laser amplitudes. In Fig || the lasers are 
weak and we see two narrow peaks from the two transi- 
tions. In Fig. ^ the lasers are ten times stronger com- 
pared to the first case. The energy levels are split and 
the spectrum has got an eight peak structure, which can 
be understood as follows: Because of the strong lasers 
all three energy levels are split into sublevels. At low 
intensity we have observed two spectral components, at 
higher driving fields each one is seen to split into four 
subcomponents. The frequencies are different for all four 
transitions. This corresponds exactly to the number of 
degrees of freedom in a normalized 3x3 density matrix. 
The two time averages have an oscillating structure when 
the time separation r is large. In practical calculations, 
infinities in equation (^) have been replaced by large time 
values. As a result we get asymmetries in the spectra and 
also finite lincwidths. 



stationary spectrum 




frequency 



FIG. 3. The stationary spectrum for our three level atom 
when lasers are weak. The energy levels and decay con- 
stants are the same as in Fig. 0. The laser amplitudes are 
0,1=0,2 = 0.2. 



C. The physical spectrum 

One weakness of the WK-spectrum is that it is only 
defined if the system is in a stationary state. This is 
clearly a limitation. For example, if instead of having 
a driving laser with constant intensity, we took a laser 
whose intensity changes in time, we obtain fluorescence 
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g-{rf-iujf}{t-ti)^-{rf+iujf)(t-t2)(y*^-^^^Y^^^-j^ 



stationary spectrum 




frequency 



FIG. 4. The stationary spectrum for our three level atom 
when the lasers are strong. The energy levels and decay con- 
stants are the same as in Fig. 0. The laser amplitudes are 
0,1 = 0.2= 2.0. 



light with varying intensity and spectrum. The fluores- 
cence light displays a spectrum but we cannot calculate 
it using the WK-definition. There are several suggested 
definitions of a spectrum which would be physical also 
for nonstationary systems In the stationary limit 

all these spectra give the WK-spectrum. 

A generalization which takes into account also the 
measurement scheme was proposed by Eberly and 
Wodkiewicz What they call the physical spectrum, 
is defined as 

oo oo 

SpHYsit^ujf.Tf) ^ J dh J dt2H*{t-h,ujf,Tfy 



— oo — oo 



H{t-h,LOf,Tf){V*{h)Vih)), (10) 

where H{t, ujf,Tf) is a filter function characteristic of the 
filter (see Fig.0) 



oc 



*duj. 



(11) 



H{uj,ujf,Tf) is a function of lo, the parameter Uf is the 
mean value of the passband of the filter and is its 
width. The filter function of a Fabry-Pcrot filter, which 
is used in our calculations, has the form 

H{t,u;f,Tf) = Qit)Tfexp{~{Tf + tLOf)t), (12) 

which gives for physical spectrum the expression 



z z 

SpHYs{t,i^f,Tf)^Tj J d.ti J 



Ff / dri I dT2e 



{V*lt-T,)Vlt-T2)). 



(13) 



dU 



From the definitions above, we see that as well as for 
the Wiener-Khintchine spectrum, two time averages are 
needed. The difference compared with Swk is that now 
we need only correlation functions whose times are less 
than the time at which we calculate the spectrum. 



1. Calculation of the physical spectrum using QRT 

In order to calculate the physical spectrum at time 
T, we need to calculate the correlation functions 
{dHti)d{t2)), < ti,t2 < T. For numerical calcula- 
tions we have to discretize the time ii, t2 = n ■ At, n = 
0...iV, Ai = There are now N time instants in be- 
tween t = and t — T. The algorithm is the following: 

1. At t — 0, form Cij{t2 = 0,ti = 0) using the initial 
density matrix g{0) and use the algorithm presented in 
the last section to calculate Cy {t2, ti) for the time values 
t2 = 0, < ti < T. 

2. Evolve g{t) to the time At, form Cij{At, ti) and use 
QRT to evolve it to ti = T. This gives the correlations 
at t2 = At, At < ti < T. 

3. Evolve g(t) to the next time step, form with a 
new t2 and evolve it to ti — T. Repeat this step until 

t2=T. 

4. The algorithm above gives correlations for time 
pairs t2 — n ■ At, t2 < ti < T, ie. t2 < ti. In or- 
der to obtain results for time values in the reverse order, 
take the complex conjugate C(ti,t2)* = C(t2,ti). 

Once we have calculated the correlation functions on a 
grid dense enough, the physical spectrum is a sum over 
the grid according to equation (|l^) . 

The detection operator used is the same as in 
the stationary case (^. Now the expectation values 
(0^(ti)0(t2)) are calculated using the algorithm for a 
nonstationary spectrum. Figures |^,^ and |^ show time 
dependent physical spectra calculated by the method de- 
scribed above. The parameters are the same as in Fig. 
^, so we expect the steady state spectrum to show eight 
peaks. The initial state of the system is g(0) — |2)(2|. 
For the time dependent spectrum also the initial state 
has some influence. Spectra with different initial states 
are different for short times. 

From Fig ^ figure it is seen that when time is small 
the spectrum has two broad peaks. At t = 4.0 a little 
more structure begins to appear, and at t = 8.0, eight 
peaks can be seen. When time increases the spectrum 
approaches the steady state spectrum, cf. Fig. ^. The 
broad spectrum features at small times derive from the 
uncertainty principle, which does not allow us to see too 
detailed structures of the spectrum, or more precisely. 
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the spectrum does not display any detailed structure at 
small time values. 



Physical spectrum 
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2. The theory of cascaded open system 

Next we show how the theory of cascaded open sys- 
tems |l^Jll|] can be used to calculate the time dependent 
spectrum. By this we mean a system (B) which is driven 
by radiation coming from another quantum system (A). 
Radiation from the system A, which may be for example 
fluorescence light, is mediated through a reservoir to a 
system B. The coupling is uni-directional, the radiation 
from system B does not go to the system A. The system 
configuration is shown in Fig. ^j. In the following, we 
derive a master equation for such a system, following the 
presentation in the paper by Carmichael [0 . 

The initial Hamiltonian is the following 



H = Ha + He 



H 



AR 



H 



BR, 



(14) 



where Ha and Hb are the Hamiltonians for systems A 
and B respectively. Hj^ is the reservoir Hamiltonian 
which mediates radiation from A to B. Har and Hbr 
are coupling Hamiltonians to the reservoir 



HAR = iVlA{aE\0)~h.c) 
HBR^iV^ibEHl)-h.c), 



(15) 



where a and b are annihilation operator for system A 
and B respectively. The coupling to the reservoir is in 
two different spatial locations for the systems A and B. 
We next proceed to eliminate this complication. Using 
the unitary transformation 



FIG. 6. The time dependent physical spectrum calculated 
using QRT at the times t = 4.0 and t = 8.0. The parameters 
the same as in Fig. ^. 
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FIG. 7. The time dependent physical spectrum calculated 
usin^QRT at time t = 16.0. The parameters the same as in 
Fig. 



Ua{t) = exp[i{HA + Hr + Har)t], 



(16) 
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we obtain an equation for the transformed density matrix 
(J with Hamiltonians which are now in the standard form 



H' — Hq + Hn + H; 



where 



Hs = Ha 
HsR i{^/jAa 



SR, 



Hb + -j^^/lAlB{a)b - h.c) 



7^S)(£;t(o)-h.c) 



(17) 



(18) 



The master equation can now be derived in the conven- 
tional manner by adiabatic ehmination of the reservoir. 
With these Hamiltonians the derivation gives 



Cg'C^ 



^g'C^C, (19) 



where the Lindblad operator C is 



(20) 



Next we give a short description how the theory de- 
scribed above can be used for spectrum measurements. 
As our detector system we choose a two level atom which 
is in the ground state. When we start to drive the atom 
using radiation with very low intensity, the excited state 
of the atom appears with a small probability. The proba- 
bility of the excited state is a function of intensity and the 
spectrum of the radiation. If the atom has a very narrow 
linewidth, then only that part of the radiation which is 
in resonance can excite it. Because excitation is linearly 
proportional to the low intensity, the probability of the 
excited state is proportional to the the intensity of the 
incoming radiation at the atomic resonance frequency. 





FIG. 8. The schematic layout of cascaded open systems. 
Fluorescence radiation from atom A goes to atom B, but the 
fluorescence radiation from B does not go back to the system 
A 



used to analyze the radiation emitted from the system 
A. The atom B has a small decay constant Tg which 
corresponds to F^? in the definition of the physical spec- 
trum (|l0|). The system A is driven by a laser. Only 
a very small fraction of the fluorecent radiation coming 
from system A is guided to atom B ic. the constant p is 
very small. The majority of the radiation (1 — p)Ta goes 
into other directions and misses the detector. Division 
of the radiation is needed, because otherwise the incom- 
ing radiation would saturate the atom B and excitation 
would not be linearly proportional to the intensity any 
more. Figure |l^ shows the probabilities of the excited 
state of an assembly of analyzer atoms as functions of 
time for the system described in detail later in this paper. 
Different curves correspond to different values of the res- 
onance frequency. The spectrum is obtained by plotting 
the excitation probability at any given time with lob as a 
parameter. This method does not need to evaluate mul- 
titime averages as in all earlier definitions. The spectrum 
is obtained from one time averages of another quantum 
system, which is coupled to the system we are interested 
in. This method of spectrum calculation is quite close to 
the method presented in the book by M.Sargent III et.al. 

@- 

The system studied was described in section II. The 
time dependent spectrum of this three level atom has 
been calculated using a two level atom as a spectrum 
analyzer. The master equation for the three level atom 
and an analyzer atom according to the theory of cascaded 
open systems is the following. 



dg 



-i[HA + HB+Hc,g]- 

4 



. . 1 . . 1 „ 
'^{UgLl - -L\Ug- -gL]U), 



(21) 



1=1 



where 



Ha - + (|l)(2|e'(-^-i)* +h.c) + 

i=l 

(|2)(3|e'('^^"'^^)*-Hh.c) 



Hb — -lobctz 



(22) 



He = - V7lpr^(|2>(l|<T- - h.c) + 



y72prB(|3)(l|a_ - h.c) 



III. CALCULATION OF THE TIME DEPENDENT 
SPECTRUM FOR THE THREE LEVEL ATOM 

A. The master equation for spectrum analyzer 

The configuration for a spectrum analyzer is shown 
in Fig. ^. The system B is a two level atom which is 



and the Lindblad operators are 



Li = ^Tai{1-p) |1)(2| 
L2 = V^A2{1-P) |1)(3| 
^3 - VtI^\1){2\ + VO.SFs 



(23) 
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In addition to the atomic Hamihonians, there is an ad- 
ditional term in the Hamiltonian arising from the deriva- 
tion which gives a small energy shift. The first two decay 
operators describe decays of the three level atom which 
take account of the majority of the radiation. They cor- 
respond to the (1 — p)Ta term in Fig. ||. The last two 
terms describe that part of the radiation which goes to 
the analyzer atom, the pTa term in Fig. ^ They also 
include spontaneous emission from atom B. 

It can be seen that that if p = the atoms are not 
coupled. Then, if the initial state is factorizable then 
the two atoms evolve separately. Atom A has two decay 
terms, which describe the transitions |1)(2| and |1)(3|. 
The Analyzer atom has one decay term with the decay 
constant F b ■ This separation is natural because if p = 
there is no radiation coming from atom A to atom B. 
If p is nonzero it is possible to trace atom B out of the 
equations and we get the well known master equation for 
atom A. This means that system B does not affect the 
time evolution of system A at all. 




pi"* 




FIG. 9. This is the version of the cascaded open system 
we use as a time dependent spectral analyzer. A small portion 
p of radiation from atom A goes to atom B. The majority of 
the radiation goes into other directions (1 — p). The system 
B decays with decay constant Fb- 



B. Simulations studied 

The first case to be studied is our three level system 
with weak lasers. Figure ^ shows the time evolution of 
the diagonal elements of the density matrix. The initial 
state is g{0) = |2)(2|. At first the population of level 
two decreases and starts to Rabi oscillate, at large times 
it approaches a constant and reaches its steady state. 
The populations of levels one and three increase in the 
beginning and for large times reach steady state. The 
behaviour of the off-diagonal elements is quite similar 
(Fig. |l^) , first there are oscillations which smoothen out 
at large times. In Fig. |l^ the population of the excited 
level is shown for all 128 analyzer atoms used in our simu- 
lations. The different curves belong to different analyzer 
atoms with different frequencies. They are clearly excited 



differently and this means that the atoms can really rec- 
ognize the spectral structure of the incoming radiation. 
The simulations were done using the Monte-Carlo wave- 
function technique ||g|,[l3|-|l5|] . The curves are not very 
smooth so bigger ensembles would be needed. The nu- 
merical calculations were, however, very demanding and 
128 different target atoms and the size of their ensembles 
were chosen for practical reasons. 

Spectra read from the analyzer atoms with weak lasers 
are shown in Fig. 13|, |lj and ^ The parameters are 
the same as in Fig The parameter p in the Lindblad 
operators ( |2^ ) is chosen to be p = 0.005. The decay 
constant for analyzer atoms is Tb = 0.001. The initial 
state is ^(0) = |2)(2|. 

In figure [l^ the spectrum is broad. At t = 2.0 there 
is a broad peak centered around w = 4 and as time in- 
creases the peak becomes narrower which can be seen at 
t = 8.0 (Fig. p^ ). The appearance of a peak around 
uj — 4: is understandable, because initially level two was 
populated. The spectrum reveals that the atom decays 
from level two to the ground level. As time increases, 
another peak begins to appear at w = 8.0. Its relative 
intensity increases and in steady state it is bigger than 
the oj = 4.0 peak. At long times, t = 200, the spectrum 
reaches its steady state Fig. ^ cf. Fig ^. 

The spectra with stronger lasers are shown in figures 
P^ , |l7| and |l^. At small times, the spectrum is broad and 
the more precise structures appears when time increases. 
In the last figure the time is very large t = 200 and 
spectrum has reached its steady state Fig. |[ 



Excitation of analyzer atoms 




FIG. 10. Excitation of all 128 analyzer two-level atoms 
as functions on time. For most of the atoms, the excitations 
is so small that it cannot be seen. The three level system 
parameters are the same as in Fig. 0. The laser amplitudes 
are Qi = Q2 = 2.0, the parameter p = 0.005, and the decay 
constant of the analyzer atoms Fb = 0.001. Each curve in 
this figure is obtained from an ensemble of identical atoms 
with 300 members. 
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Diagonal elements 




10 20 30 40 50 60 70 80 90 100 

time 



FIG. 11. Time evolution of diagonal elements. The energy 
levels and decay constants are the same as in Fig. 0. The laser 
amplitudes are Qi — ^2 = 2.0. Dashed line: gn, solid line: 
Q22, dashdot line: £133. 
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FIG. 13. The time dependent spectrum read from the 
analyzer atoms at times t = 1.0 and t = 2.0. The parameters 
are the same as in Fig. ^ 




16 



FIG. 12. Time evolution of absolute values of off-diagonal FIG. 14. The time dependent spectrum read from the 

elements. The parameters are the same as in Fig. |ll|. Dashed analyzer atoms at times t = 4.0 and t = 8.0. The parameters 
line: 1^12!, solid line: \qi3\, dashdot line: 1^231. ^''^ ^^"^^ ^ i° Fig. § 
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Spectrum read from analyzer atom 



Spectrum read from analyzer atoms 
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FIG. 15. The time dependent spectrum read from the 
analyzer atoms at time t = 200.0. The parameters are the 
same as in Fig. ^ 



FIG. 17. The time dependent spectrum read from the 
analyzer atoms at times t — 4.0 and t — 8.0. The parameters 
are the same as in Fig. ^. Compare this with Fig. g calculated 
with the same parameters. 



Spectrum read from analyzer atom 




Spectrum read from analyzer atoms 



10 12 14 16 
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FIG. 16. The time dependent spectrum read from the 
analyzer atoms at times t — 1.0 and t — 2.0. The parameters 
are the same as in Fig. ^. Compare this with Fig. Q calculated 
with the same parameters. 



FIG. 18. The time dependent spectrum read from the 
analyzer atoms at time t — 200.0. The parameters are the 
same as in Fig. ^. Compare this with Fig. ^ calculated with 
the same parameters. 
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C. Comparison between methods 

The spectrum calculated using analyzer atoms and the 
definition of the physical spectrum seem to give very sim- 
ilar results. However, the calculation techniques are to- 
tally diff'erent and from the mathematics it is not at all 
obvious that these two spectra should agree. For the 
physical spectrum, two time averages are needed and the 
spectrum is a kind of Fourier-transform of them. In this 
respect the physical spectrum is closer to the Wiener- 
Khintchine spectrum. The method with analyzer atoms 
does not need any multitime-averages. The spectrum is 
read from a quantum system (the two state atom) which 
acts as a measurement device. 

It is also interesting to compare how the uncertainty 
principle restricts spectral measurements. In all our sim- 
ulations, analyzer atoms of very small linewidth were 
used. This means that the atoms respond to the changes 
of the incoming radiation very slowly, but if the interac- 
tion time is long enough the atoms give an accurate spec- 
trum. If we used atoms with a larger linewidth, the atoms 
would indicate the state of the system more quickly, but 
because of their large linewidth they would not give a 
very accurate spectrum. In the case of a the physical 
spectrum, the uncertainty relation is related to the prop- 
erties of the Fourier transform. 

The method of analyzer atoms is based on the sim- 
ulation of a master equation. Recently many paral- 
lel stochastic algorithms have been presented [|ll|,|l^-^5j . 
There are also freely available programs which are devel- 
oped for simulations of these kinds of master equations 
| fi6[ . It is also possible to parallelize the algorithm by 
computing the evolution of the different analyzer atoms 
on different processors. 

IV. CONCLUSION 

We have shown a fundamentally new method to cal- 
culate a spectrum. Instead of calculating two-time aver- 
ages using the Quantum Regression Theorem we couple 
the system to two state atoms which are used as a mea- 
surement device. The spectrum is read from the exci- 
tation probability ie. from one-time averages of the two 
state atoms. This method allows us to calculate both the 
time dependent and stationary spectrum. Comparisons 
between the spectra calculated show that using this new 
method we find similar results as using two-time averages 
and a Fourier-transform. 

Because of the frequency-time uncertainty, the spec- 
tral resolution is poor for initial times. However, when 
the bare level structure is considerebly modified by the 
presence of strong laser fields like in Fig. ^ it may not 
be possible to surmise the low intensity structure. Fig. 
^, which is characteristic of the system of interest. The 
the short time spectra, Figs. ^ and |l^ may indicate the 
number of transitions involved even if their positions (4 



and 8 in our case) are very difficult to locate. When 
the measurement goes on, more data are collected and 
the spectral information approaches the steady state as 
closely as we desire for a meaningful measurement. 
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